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1 Introduction 

Wavelets were introduced in the beginning of the 1980s synthesizing ideas orig- 
inated in engineering (subband coding in signal processing, pyramidal and mul- 
tipole algorithms in image processing), physics (coherent states, renormaliza- 
tion group), numerical analysis (spline approximation) and pure mathematics 
(Littlewood-Paley and Calderon-Zygmund theories), and provide us with a par- 
ticularly simple and powerful tool in mathematical analysis with a great variety 
of applications [3 [S [H [13] . 

Two unitary operators on L^(M) play an important role in wavelet theory: 
dilation by r > 0, Dr, defined by [Drf]{x) := r^/^ f{^x), and translation by t, 
Tt, defined by [Ttf]{x) := f{x — t). Wavelet bases are built by application of 
translations and dilations to an appropriate function. Since the main stream 
of development leading up to wavelets departs from Fourier analysis, Fourier 
transforms and series are involved in most of the techniques in wavelet theory. A 
clear reason for this is the fact that the complex exponentials are eigenfunctions 
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of the derivative operator and then give us an spectral representation of the 
unitary one-parameter group generated by it, i.e., the group of translations 
{Tt-.te R} on L2(r) [J. 

Paying attention to the discrete theory of wavelets, the values of r and t are 
fixed (usually r = 2 and t = 1) and one works with the entire powers of the 
operators D := D2 and T := Ti. The spectra of D and T coincide with the unit 
circle 9D of the complex plane C and have constant denumerable multiplicity. 
Thus, the functional spectral representations of both operators live in direct 
integrals L^(i9D,7i) of valued functions over 9D, where H is an auxiliary 
separable Hilbert space [TJ [71 [T7] . Direct integrals have been already considered 
in abstract wavelet theory -see [13] and references therein-. Here all these facts 
are proved in Section [2l Explicit spectral representations of D and T are given 
in Propositions [1] and [2l They are built on the basis of orthonormal bases of 
L^(M) with appropriate structure and are useful for computational purposes. In 
practice, any orthonormal basis of i^(K), after adapting it, can be used to build 
such spectral representations. To illustrate this the details for the exponential 
and Haar bases are included in an Appendix at the end of the paper. Moreover, 
an additional spectral representation for T is defined in Proposition [4] on the 
basis of the usual Fourier transform. Remark [5] and Proposition [1^] evidence 
that classical results in wavelet theory join our approach under this spectral 
model. 

The main purpose of this work is to characterize orthonormal wavelets and 
multiresolution analysis (MRA) of L^(R) by means of rigid (operator- valued) 
functions defined on the functional spectral spaces of D and T. A rigid function 
on L^(c)D, Ti.) is a function A from into the space >C(7i) of bounded operators 
on Ti. such that A{u!) is a partial isometry with the same initial subspace for 
almost every (a.e.) oj G 9D. For orthonormal wavelets this is done in Section 
[41 fTheorem [TU]) and for MRA in Section l5l fTheorems ITSl and l20l) . The key to 
obtain these characterizations is the close connection between wandering and 
invariant subspaces of L'^{dV),7{) and rigid and range functions defined on it 
[H [51 [m [TS] . The fundamental concepts and results of this theory are collected 
in Section [31 

Section 15.11 is devoted to get the pair of discrete quadrature mirror filters 
associated with an MRA from the two-scale relations in this context. Remark 
[26l makes clear these filters act just on the generic spectral models given in 
Proposition [21 

Theorem [121 Example [TU Corollary [191 and Proposition [281 highlight the 
usefulness for computational purposes of the spectral models considered in this 
work. Finite sets of non-zero coordinates imply local character in "time" and 
"frequency" , the meaning of these concepts depending on the bases considered. 
At the same time, suitable choices of the bases lead to sparse sets of conditions 
(compare, for example, the matrices (a*'//") for the exponential and Haar bases 
in Appendix). A detailed analysis of this set of conditions shall be done else- 
where. 
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2 Spectral models of dilations and translations 



Let D denote the open unit disc of the complex plane C and dD its boundary: 

D := {A e C : |A| < 1}, dB := {uj e C : ^ 1} . 

In interpret measurability in the sense of Borel and consider the normalized 
Lebesgue measure du!/{2iT). Given a separable Hilbert space H, let L'^{dIi]Ti) 
denote the set of all measurable functions v : 3D H such that 

J da 

(modulo sets of measure zero); measurability here can be interpreted either 
strongly or weakly, which amounts to the same due to the separability of Ti.. 
The functions in L^(9D;7i) constitute a Hilbert space with pointwise definition 
of linear operations and inner product given by 

(u,v):=/ (uH,vH)^|^, (u,veL2(5D;7^)). 

J do 

The space L'^{dB;7i) is a particular case of direct integral of Hilbert spaces. 
The theory of direct integrals is originally due to von Neumann [TT] and is in 
the basis of the functional spectral models for operators [U Chapter 7]. In the 
case of a constant field of Hilbert spaces, Hui = H for all uj g 9D, the direct 
integral J^Ti-cu |^ is just the space L^{dIS>;H). Elementary properties of vector 
and operator valued functions can be found in [71 Chapter HI] ; see also [T and 
[g Lecture VI]. 

Now, let r > 1 and consider the dilation operator Dr : L^(M) L^(R) 
defined by 

[Drf]{x):=r'/yirx), (feL^m). 

Dr is unitary on and then its spectrum is included in dB. In order to give 

the spectral resolution of Dr, we must find a direct integral of Hilbert spaces 
/J^ '^i^ '^/^('^) and a unitary transformation Q from L^{R) onto J^Hi^ diJ.{uj) 
such that 

[gDrf]{uj)^uj[gf]{uj), fora.ewe9Dand/eL2(M). (1) 

In other words, QDrQ^^ is the operator "multiplication by lj" in the direct 
integral: 

[gA-g-if](tj) = • f(w) , for a.c e a© and f e / ^ . 

Jaa 27r 

For it, let {K'^±'-{x)}j^j be an orthonormal basis (ONB) of L^[±l, ±r), where JJ 
is a denumerable set of indices (usually, J = N or J = Z) , so that the set 

{K^{x) :^r™/2<)(r'"x)},ej 
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is an ONB of L'^[±r-'^, foj. gvery to G Z. For each / G L'^{^) put 



7(m) 



/(a;) (x) dx , (j G JJ, m G Z) , 



Let denote the Hilbert space of sequences of complex numbers (cj)jgj[ such 
that X^jGJ I'^il^ < ^"^"^ {^''jls^i ^ fixed ONB of ^^(JJ)©P(J), where 

© denotes orthogonal sum. 



Proposition 1 The operator Q defined by 



g. L\i 



dcu 



[/'(J) © — = L^{dB; /2(JI) © P(J)) 



(m) 



(2) 



determines a functional spectral model for the dilation operator Dr, i.e., Q is 
unitary and satisfies (QP. 

Proof: Q is unitary since {k''^\x)} s=±,j^],,mi^z is an ONB of L^(R) and 
f|P = l/ij^ for every / G L2(r). Moreover, satisfies (P 



'{rn) 



because [Drf]^ j — fl"j for all s — ±, j e I, m e Z. □ 
In a similar way, for t > and the translation operator Tt on (M) defined 



by 



[TJ]{x):=f{x-t), if EL' 



given an ONB {Lf\x)}i(zi of L^[0,t), with I a denumerable set of indices, the 
translated set 



{4-\x) :=Lr'{x-tn)},^ 
is an ONB of L'^[tn, t{n + 1)) for every n £ Z. For each / G L^(IR) write 



(0)/ 



J i 



t{n+l) 



f{x)L^^\x)dx, {iel,neZ) 



and let {wij^^j be a fixed ONB of ^^(I). 
Proposition 2 The operator J- given by 



' l\l)'^=L^{dW,l\l)) 



(3) 



n t{n) 
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determines afunctional spectral model for the translation operator Tt, i.e., T is 
unitary and satisfies 

[TTtf] {uj)^uj [Tf] (tj) , for a. ecu € dB and f e (R) . (4) 

Proof: The operator T is unitary since {-Z^j-"''(a;)}iei,riez is an ONB of L'^(M.) 
and ll/ll = ||f||2 = J2i,n \ft^\^ for every / e L^{R). T satisfies © because 

\Ttf\, = /i"" ' for alH e I, n e Z. □ 

Since {i4:^™^(a;)}s=±jeJ,mez and {L|"''(x)}ieL„ez are ONB of given 
/ e one has (in i^-sense) 

/=E/i?^S^ and / = (5) 

s,j,m i,n 

The change of representation between the two models is governed by the matrix 
so that 

r^(m) _ sj,m ^ (n) r (") _ \ ' s,],m r^{m) 

i,n s,j,m 

and 

In what follows we shall focus attention on the values r — 2 and t — 1, i.e., 
the operators 

:= D2 and T := Ti . 

For these values, the elements of the matrix (a*';^'™) for the exponential and 
Haar bases are given in the Appendix. 

Fixed the orthonormal bases {-fi^±™^ (a;)}jej,mGZ and {i|"^(a;)}igi_„gz of L^(M) 
giving rise to the respective spectral models ^ and ^ oi D and T, for / G 
L^(]R) we shall write 

Tf = i = {ft^}, g/ = f = 

Lemma 3 Let f G L^(M) anrf, for p, q G Z, 

g := DPT^ f, h := T'^D^' f . 

If^f = {t^}, Tg = {gf)}, Th = (hf^), Qf = {/ij)}, = {gg)} 
Qh=\h^^^), then 

al"^ = E <r E ""n+T" f^"^ > (/ G I, G z) , (7) 
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= E E <nZ, fS', (r = ±, / e J, fc G Z) , (9) 

2,n s,j.7n 



h^S = E "^.'n' E /iJ^ ' (r = ±, ; e J, fc e Z) . (10) 

Proof: The result is a straightforward consequence of the change of repre- 
sentation formulas (O and the fact that QDQ^^ is the operator "multiplication 
by w" on QL'^{E) and TDJ^~^ is the operator "multiplication by lo" on TL/^iM). 
□ 

Finally, let / denote the usual Fourier transform of a function f ^ (M) : 
/(y):- / f{x)e~^^'^^ydx, (yeM). (11) 



Proposition 4 Let {wfel^.g^ a fixed ONB o/P(Z). T/ie operator T^, defined 
by 



where, if UJ = e^''*^, 



(12) 



/fc(c^) = /fc(e2"'') := /(^ + k) , /or a.e. G [0, 1) and k £ Z , 
determines a functional model for the translation operator T . 
Proof: The operator is unitary since, for / G L^(M), 

ii/ir - ii/ip - E i/> +k)?de^Y.I \M 

Moreover, if w = e^'^'^, then 



ZTT 



[T/]fc(^) - [T/](e + fc) = e-2.^(e+fc) j(0 + fc) = . /fe(^) , 
so that [J^,T/](w) = w • [J^,/](w), for a.e. w G a© and / G L^{R). □ 
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3 Wandering and invariant subspaces. Rigid and 
range functions. 

Let H he a separable Hilbert space and denote by C{Ti,) the space of bounded 
linear operators on H. Consider the functional space L'^{d3;TC) defined in 
Section [2] and the subspace £ of H) consisting of all constant functions, 

i.e., the functions v : dU) TL such that there exists a vector v E H with 
v(ti;) = V for a.e. lo G dUi. From now on the operator "multiplication by oj" on 
L'^idB; n) shall be denoted by M, that is, 

[Mv](a;) :=wv(w), e L^{dB ■,?{), uj e dB) . (13) 

It is easy to verify that M is unitary and M^^{— M*) is defined by [M*v]{lu) = 

LU* ■ V(w). 

A subspace of a Hilbert space is called a wandering subspace for an operator 
U if it is orthogonal to all its images under the (positive) powers of U. For 
isometrics, wandering subspaces behave better than usual: if U is an isometry 
and if 371 is a wandering subspace for U, then [/™97l _L [/"9Jl whenever m and 
n are distinct non-negative integers. If U is unitary, even more is true: in that 
case U^VJl _L U^VJl whenever m and n are any two distinct integers (possibly 
negative) . 

A weakly measurable^ operator- valued function 

A:dB^ C{n) : lu ^ A{u;) 

is called a rigid operator function if A(uj) is for a.e. lu G (9D a partial isometrjf^ 
on Ti. with the same initial space. 

According to Halmos [4j Lemma 5], wandering subspaces for M and rigid 
functions are related as follows: 

Lemma 5 [Halmos] A subspace 9Jl of L'^{dO] TL) is a wandering subspace for M 
if and only if there exists a rigid function A such that DJl ~ A€. The subspace 
971 uniquely determines A to within a constant partially isometric factor on the 
right. 

On the other hand, a closed subspace 971 of i^(9D; H) is called invariant if 
Mv G 971 for every v G 971. 971 is called doubly invariant if Mv and M~^v belong 
to 971 for each v G 971. 971 is called simply invariant if it is invariant but not 
doubly invariant. 

A range function J = J{u!) is a function on dH taking values in the family 
of closed subspaces of Ti. J is said measurable if the orthogonal projection P{lo) 
on J{uj) is weakly measurable. Range functions which are equal a.e. on (9D 

-'That A is weakly measurable means that the scalar product {A{u))h, g)ti is a Borel mea- 
surable scalar function on 9D for each h,g ^Ti.. 

^An operator B £ CCH) is a partial isometry if there is a (closed) subspace of W such 
that ||-Bn|| = ||«|| for « G 9Jt and Bv = for v £ SXTt^. In such case 9Jl is called the initial 
space of B. 
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are identified. For each measurable range function J, dJlj denotes the set of ah 
functions v in L'^{dIS);'H) such that v{lo) lies in J{w) a.e.. The correspondence 
between J and Tlj is one-to-one, under the convention that range functions are 
identified if they are equal a.e.. 

The following result is implicitly contained in the work of Lax [TTI and 
explicitly proved by Srinivasan [TS] -see also Theorem 8]-. 

Lemma 6 [Lax-Srinivasan] The doubly invariant subspaces of L'^{dD;Ti,) are 
the subspaces where J is a measurable range function. 

In what follows by the range of a doubly invariant subspace 37lj we mean 
that range function J. This definition can be extended to an arbitrary set of 
functions: the range of a set of vector functions is the range of the smallest 
doubly invariant subspace containing all the functions. 

To determine the simply invariant subspaces of L^(9D; TC) one needs to in- 
troduce the Hardy classes. We denote by i?^(D; H) the Hardy class of functions 

oc 
k=0 

with values in H, holomorphic on D, and such that ^ J^j, ||u(rtj)||^ duj, (0 < 
r < 1), has a bound independent of r or, cquivalently, such that ^ < oo. 

For each function u e i/^(B;7Y) the non-tangential limit in strong sense 

oo 

s-limu(A) = '^^uj'^ak =: u(a;) 

k=a 

exist for almost all uj € 9D. The functions u(A) and u(a;) determine each other 
(they are connected by Poisson formula), so that we can identify if^(D;7Y) 
with a subspace of L2(9D;H), say H+{d]D);n), thus providing H'^{]D);n) with 
the Hilbert space structure of 77+((9D;7i) and embedding it in L^(9D;7i) as a 
subspacell The complementary space H-{dB);n)j_= L^{d]0);H)\H+{dD;n) is 
associated with the conjugate Hardy class 7?^(C\I]); 7i) in a similar way. 

The general form of the simply invariant subspaces of L^(i9D; H) is given by 
Helson [5l Theorem 9] on the basis of the work of Halmos [¥, Theorems 3 and 
4]. 



Fixing some orthonormal basis {«! , 1*2 1 ■ ■ ■} for H, each element f of L^(c?ID); ?-{) may be 
expressed in terms of its relative coordinate functions fj, 

j 

where the functions fj are in L'^{dn) := L^{dIi;C) and [[flp = J2j ll/jlP- A function f in 
L'^{dO;'H) belongs to _ff+(9ID); H) iff its coordinate functions fj are all in the Hardy space 
H+{dD) := H+{dn;C). 
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Lemma 7 [Halmos-Helson] Each simply invariant subspace 971 of L'^{dI};TC) 
has the form 

AH+{dI]>;H)(SmK , (14) 

where K is a measurable range function and A is a rigid function with range J 
orthogonal to K almost everywhere. The doubly invariant part 971^ of 9Jl is just 
Hn^zM^dJl. The subspace 9H uniquely determines the rigid operator function A 
to within a constant partially isometric factor on the right. 

4 Orthonormal wavelets on M 

An orthonormal wavelet on M is a function ■0 G i'^(]R) such that {ipj^k '. j,k ^ Z} 
is an orthonormal basis of L^(M), where 

ipj.kix) [D^T^ij]{x) = 2^/^(2^ a; - k) , (j, /c e Z) . 

Remark 8 It is a weU known result -see, for example, Lemma 2.2.4]- that 
for if) S -L^(IR.), {'(/'(• — k) : fc G Z} is an orthonormal system if and only if 

Sfcez + = 1 foi' ^-^^ 6* e E, where tj) is the Fourier transform of tp 
defined in pT]) . This result fits in with the concept of rigid function with one- 
dimensional initial subspace on the spectral model for T given by HH). Indeed, 
{■)/'(• — /c) : /c G Z} is an orthonormal system in L^(R) if and only if the closed 
subspace spanned by xjj^ := T^ip is a wandering subspace for the unitary opera- 
tor M — T^TT^^ defined in jr^L^(M) by IT^ . By Lemma [51 this is equivalent 
to the existence of a rigid function A on ^*L^(R) with one-dimensional initial 
subspace and such that tp^{uj) = A{u;) u, with u a normalized vector belonging 
to the initial subspace of A{u!), the same one-dimensional subspace of /^(Z) for 
a.e. u! e 9D. Since A is a partial isometry for a.e. w = e^'^*^ e 915, 

1 = Iklli^(z) = II^Hwlli^(z) - = 

= ^|^Afc(e2"'^)P = ^|V3(e + fc)|2, for a.e. 0€R. 

fcez fcez 

From now on in this Section we shall pay attention to the respective spectral 
models for D and T given in Propositions [1] and [2l £ shall denote the subspace 
of C/L^(K) = L'^{dB; P(J) © P(J)) of all constant functions and £ the subspace 
of constant functions of TL'^{R) = L'^{dl]i; that is, 

£ g(L2(-r, -1] © ^^[1, r)) , £ := ^'^(^^[0, 1)) . 
The expression (•) shall denote the closed subspace expanded by the argument 

(•) := span{-}. 

We shall use the discrete Dirac delta function 6k defined by 

„ _ / 1, if fc = 0, 

{ 0, if k e Z\{0}. 
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The next result is straightforward. 

Lemma 9 Let tp be an orthonormal wavelet of L'^iM). Then: 

(a) {ijj) is a wandering subspace for D and T on L^(M). Equivalently, Q{ip) 
is a wandering subspace for QDQ^^ on OL'^(R.) and J-'{ip) is a wandering 
subspace for J-TJ-~^ on TL'^iM.) 

(b) {T'^iIj : k ^ Z) is a wandering subspace for D on L^(R). Equivalently, 
Q{T^^] : fc e Z) is a wandering subspace for QDQ^^ on t/L^(R). 

(c) Q{T''ip : fc e Z) is a doubly invariant subspace of TL'^(K). 

Now, the resuhs of Section [3] aUow us to characterize any orthonormal 
wavelet of L'^(R) in terms of rigid functions on gL'^{R) and TL^(R). 

Theorem 10 If ip is an orthonormal wavelet o/L^(R), then there exist rigid 
functions A on J-L'^{R) and B, C on QL'^{M.) such that: 

(i) the initial subspaces of A, B are one- dimensional, and C{uj) is unitary for 
a.e. uj £ 9D; 

(ii) T-\uj''Ai: fc e Z) = g-^{C€); 

(iii) {i}) = F-^{At) = g-^{Bl). 

The wavelet '0 uniquely determines A, B and C to within constant partially 
isometric factors on the right. 

Conversely, assume that there exist rigid functions A on J-L'^{R) and B, C 
onQL'^iR) satisfying (i), (ii) and (Hi') T-^{A€) = g-^{B€). Then, the triplet 
{A,B,C) has associated a unique orthonormal wavelet ip G L^(R) given by 

^■.= J^-^[Au]=g-^[Bw], (15) 

where u(w) — u and v(u;) = v for a.e. lu E 9D, being u and v normalized vectors 
in the initial spaces of A{uj) and B{lj), respectively. 

Additionally, there exists a measurable range function J = J{uj) on J-'L'^(R) 
such that J{uj) is one- dimensional for a.e. lo G dH and 

T-^Mj = T-^{uj^Ai:ke'L)= g-^{C€) . (16) 

Proof: Let V' be an orthonormal wavelet of L^(]R). According to Lemma 
[51(a) and Lemma [51 there exist rigid functions A on TL'^i^) and B on ^L^(R) 
satisfying (iii). A and B have one-dimensional initial subspaces because {tp) 
is one- dimensional. Since TTT^^ is the operator "multiplication by w" on 
^L^(]R), Lemma [91(b) and Lemma [H ensure that there exists a rigid function 
C on gL'^{9P) verifying (ii). Moreover, CiuS) is unitary for a.e. lo G because 
gog-^ is the operator "multiplication by w" on gL'^{R) and {D^T'^ip : j,k e Z} 
is a basis of L^(R), so that 

{ij^ gT^^P : J, fc G Z) = g{D^T''^ :j,keZ)^ ^^'(R) , 
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i.e., (CC) has full range -see [51 Lecture VII]-. Obviously the function tp 
uniquely determines the subspaces (ip) and (T'^ip : k £ Z). Thus, by Lemma 
[51 ijj uniquely determines the rigid functions A, B and C to within constant 
partially isometric factors on the right. 

Conversely, let {A, B, C) be a triplet of rigid functions satisfying (i), (ii) and 
(iii'), and define V e ^^(R) by ([I5|. Then, by Lemma [51 {D^T'^i) : j,k G Z} 
is an orthogonal system because A, B and C are rigid functions. One has 
||£)j7^fc^|| = I for every j,k G Z, since D, T are unitary operators, A(u)) is a 
partial isometry for a.e. w G 9D and u G is a normalized vector. The 

completeness of the orthonormal system {D^T'^ip ■ j,k G Z} follows from the 
unitarity of C{lu) for a.e. oj G 9D as before. 

Finally, Lemma [91(c) and Lemma [6l imply that there exists a range function 
J satisfying JTB]) and such that J(tj) = {[T'iI;]{lu)) = {A{lj)u) = {B{lj)v) for 
a.e. CO edB. □ 

Remark 11 The rigid functions B and C of Theorem [TUl are related as follows: 
Let 'iTli and ^2 be the respective (constant) initial spaces of B and C. There 
exists a rigid function G on C/L^(R) = © P(JI)), with initial space 

Vti and the range of G(w) included in for a.e. uj G 9ID), such that G£ C 
P(J) ® P(J)) and B = CG. See [4, Lemma 6] for details. 

For practical purposes it is advisable to deal with coordinates. 
Theorem 12 A function ijj G ^^(R) mth Q^p ^ ij^ ^ {V^lj^} 

is an orthonormal 
wavelet if and only if the following two conditions are satisfied: 

(i) Orthonormality: For every p,q G Z, 



s,j.m i,n r.k.l 

(a) Completeness: For every finite set ¥ of indices (s,j), with s = ± and 
j G J, the maina@ 



I s,i,m "Sr^ r,k,l 7(1) \ 

\i,n r.k.l J (m,g)eZxZ 

(s,j)eF 

/las maximal range, i.e., the cardinal of¥. 

Proof: By definition, a function ip G L^(M) is an orthonormal wavelet if 
and only the system {V'p.g — D^T'^i/j : p, g G Z} is an orthonormal basis of 

*Each file of this matrix corresponds with fixed values of the pair of indices (m, 5) G Z X Z, 
whereas each of its columns corresponds with fixed values of the pair {s,j) S F, so that the 
matrix has an infinite number of files and the number of its columns is the cardinal of ¥. 
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Z/^(M). Since D and T arc unitary operators, the orthonormality of this system 
is equivalent to 

{^p.q,4>) = 5p5q, {p,qeZ). (18) 
Putting Gipp,q — 'ipp^q — {[V'p,g]ij^}i condition p8|) can be written as 

/ u;P{^,,^{i.),,p{u;))^^Sp5q, {p,q€Z). (19) 

If p 7^ 0, this means that (i/jg ^{u!),ip{uj)) , which is a complex-valued integrable 
function, is such that all its Fourier coefficients are equal to zero. If p = 0, this 
means that {^ipQ g{u!),ip{uj)) has Fourier coefficients equal to zero, except the 
constant term, which is 1. That is, ([T9|) is equivalent to 

(V'o.g(w), V'(t^)) = Sq , for a.e. uj e dB . (20) 

From 



(^o,,(^), rPico)) = E [E ^"^ [E^^M^ 




s ,j m / 

= E-^[E^i:^i^3f^ . 

Thus, in terms of coordinates, is written as 

E^i?^Or^ = '5p^,, (P,'?ez). 

Since, by ®, [^o,,]!""^^ = E,„ one gets the con- 

dition (i) of the result. 

Now, for / G (M) let us write Qf — i — { /j™^ } . The subspace 5^ consisting 
of all functions / e -L^(M) with a finite number of non-zero coordinates /^^™'' 
is dense in L^(R). Thus, the completeness of the orthonormal system {V'p.g : 
p, g e Z} in L^(M) is equivalent to the condition: if / G and (/, t/jp.q) = for 
aU p,qeZ, then / = 0. From jS]) and 

s,j.m s,j^m i.n r^k^l 

Since m and p are independent, one arrives at the condition (ii) of the result. 
□ 

Remark 13 Rearranging sums in (fT7|) and using ^ one gets 

E [^p.o]!"'' ["^o.gji"^ = (-00,9, i'pfi) = SpSq, (p, g e Z) , 
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where Tipp^q — {[V^p,?]!"^}- This gives a direct proof for the orthonormahty 
condition (i) of Theorem [121 We include the former proof as an exercise of style 
in this context. 

Example 14 For both the exponential bases of Appendix lA.ll and the Haar 
bases of Appendix IA.2I one has the following result: Let tp he a. function of 
I/^(R) with compact support included in the interval [1, 2] (any other interval of 
length one could be considered). Then, the only possible non-zero coordinates 
of ip are G J, and is an orthonormal wavelet if and only if 

i k 

and for every finite set F of indices (s, j), with s — ± and j G J, the matrix 

fE<^:;^n) 

\ fc / (m,g)eZxZ 

(s,i)GF 

has range the cardinal of F. It is easier to deal with these conditions when 
considering the Haar bases of Appendix IA.2| because the matrix (al'^J'™) is 
really sparse for them. The solutions of this set of conditions shall be studied 
elsewhere. 

5 Mult iresolut ion analysis on M 

A multiresolution analysis (MRA) of L^(R) is a sequence {V„}„gz of closed 
subspaces of L^(R) such that 

(i) {0} c • • • C V_i C Vo C Vi C ■ • • C L^iR), 

(ii) n„ezV„ = {0}, 

(iii) U„ezV„ = L^R), 

(iv) / e Vo ^ D"f e V„ for aU n e Z, 

(v) there is a function ip £ L^(R), the scaling function, whose integer trans- 
lates {T''(p : fc G Z} form an orthonormal basis of VqH 

^Sometimes (v) is replaced by the weaker condition: (v') there is a function ip G L^(R) 
such that the set {T^'ip : fc S Z} is a Riesz basis of Vq, i.e. its linear span is dense in Vq and 
there exist positive constants A and B such that ^J^fcgz'^fc — W^kez''^'^^' ~ ''^)IIl2(e) — 
^SZfcgz'^fc' ^'-'^ each {cf^} C IR such that X]fcgz'^fc < 1° fact, one can orthogonalize the 
basis and (v) and (v') are equivalent. See [12] Theorem 2.2.5], [6| Section 2.1]) or [l3l Theorem 
7.1] for details. 
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Remark 15 The properties used to define an MRA are not independent: Con- 
ditions (i), (iv) and (v) [or (v')] imply (ii) -see, for example, [6l Theorem 2.1.6]-. 
Obviously, properties (iv) and (v) imply that the set 

{ipn.k ■■= ^"TV : fc e Z} 

is an orthonormal basis of V„ for each n e Z. 

With the modest prerequisite that ip G _L^(R) and taking into consideration 
the spectral model given by in (fT^ one can simply construct MRAs: 

Proposition 16 Let ip e L^(M) n L^(R) and such that: 

(i) J2kez \'Pk{'^)\'^ = 1 for a.e. uj E dB; 

(ii) (^fc(l) = Sk- 
Then ip is a scaling function for an MRA {^n}.^^^ of L^{M.). 

Proof: According to Remark [5] condition (i) is equivalent to the fact that 
{(^(- — fc) : fc e Z} is an orthonormal system. Moreover, since (pki^^^^^) = 
(p{9 + k) for 9 £ [0, 1) and fc G Z, in terms of the usual Fourier transform p, 
condition (ii) means that <p{0) — 1 and (^(fc) = for k G Z\{0}. It is well known 
-see, for example, [TH Theorem 2.2.7]- that under these conditions the spaces 
Vn ■■= {pnM : fc G Z), n G Z, give an MRA of L'^{R). □ 

Scaling functions and MRAs of i^(R) can also be characterized in terms of 
rigid functions on the spectral models (yL^(R) and ^L^(R) given in Section [31 
The key to do it is the following simple result: 

Lemma 17 Let ip he a scaling function for an MRA {Vri}^^^ of L'^{R.). Then: 

(a) {p) is a wandering subspace for T on L^(M). Equivalently, J-{p) is a 
wandering subspace for J-TJ-^^ on J-L^{M.) 

(b) The orthogonal complement of QVq = Q{T^p : fc G Z) is a (simply) in- 
variant subspace on QL'^{M) without doubly invariant part. 

(c) TVa = T{T^ij : fc G Z) is a doubly invariant subspace of J-L'^{M.). 

Theorem 18 Letp) be a scaling function for an MRA {V„}^g^ o/L^(R). Then 
there exist rigid functions R on J-L'^(R) and S on QL'^{R) such that: 

(i) the initial subspace of R is one- dimensional and S{uj) is unitary for a.e. 
ijj G a©; 

(ii) Vo ^ T-^ijj'' R^ : k G Z) =g-i(5i7-(aD;P(J)©/2(J))),- 
(ill) {ip) =T-^{R€]. 
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The scaling function (p uniquely determines R and S to within constant partially 
isometric factors on the right. 

Conversely, assume that there exist rigid functions R on JFL-^(R) and S on 
5i^(R) satisfying (i) and 

(it') T-\uj''R€ : fc e Z) = g-^{SH-{d]D);P{^) ® P {!))). 

Then, the pair (R, S) has associated a unique scaling function ip for an MRA 
{^»}„ez ofL^m given by 

jr-i[i?w] , (21) 

where w is a (constant) normalized vector in the initial space of R. In such 
case, (a) and (Hi) are satisfied. 

Additionally, there exists a measurable range function K = K{uj) on J-'L^{M) 
such that K{uj) is one- dimensional for a.e. lu £ dB) and 

T-^mj^=T-^{uo^A€:k£l) = Vo. (22) 

Proof: Let be a scaling function for an MRA {^n} of i^(R). Ac- 
cording to Lemma 1171 (a) and Lemma [5l there exists a rigid function R on 
J^L^iM) satisfying (iii). R has one-dimensional initial subspace because (95) is 
one-dimensional. Recall that TTT^^ is the operator "multiplication by w" on 
TL'^{9:) and QDQ-^ is the operator "multipHcation by w" on QL'^{R). Then, 
Lemma 1171 (b) and Lemma [7] ensure that there exists a rigid function S on 
5L^(R) verifying (ii). Here the doubly invariant subspace OJl^ does not appear 
because n„gzV„ = {0}. Moreover, S'(w) is unitary for a.e. G 9ID> because 
and {D^T^Lp : j,k G Z} is a basis of £^(M), so that {S €.) has full range -see [H 
Lecture VII]-. Since the function ip uniquely determines the subspaces ((/?) and 
{T^(p : G Z), by Lemmas [5] and [3 R and S are determined by (p to within 
constant partially isometric factors on the right. 

Conversely, let {R,S) be a pair of rigid functions satisfying (i) and (ii'), and 
define (p G ^^(R) by Then, by Lemma El {T'^ip : fc G Z} is an orthogonal 

system because i? is a rigid function. One has ||r'''7/'|| = 1 for every fc G Z 
since T is unitary, R{u!) is a partial isometry for a.e. cu G and w G P(It) 
is a normalized vector. That U„gzV„ = L^(R) follows from the unitarity of 
S{uj) for a.e. uj G as before. That QVq has no doubly invariant part implies 
n„ezV„ = {0}. 

Finally, Lemma fTTl fc) and LemmalS] imply that there exists a range function 
K satisfying (p2)) and such that K{uj) — {[J-(p]{uj)) — {R{uj) w) for a.e. ui G 9D. 
□ 

Corollary 19 Let ip G i^(R) be a scaling function and J-ip = (p = |(^|"^}. 
Then 

^r"^ = s, . (23) 
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Proof: Condition is equivalent to ||c^(aj)|| = 1 for a.c. uj E 9©. Indeed, 



I P - E I E ^i"' r - E ( E ( E ' 



= E E ^'""'') - E E ^i"^ ^^"^) 



and the last expression is equal to 1 for a.e. uj G if and only if ([25)1 is 
satisfied. □ 

Given a MRA, denote by >V„ the orthogonal complement of V„ in V„+i, i.e. 
V„+i = W„ © V„. From (ii), 

n 

Vn+i = V„ © W„ = , (n e Z) . (24) 



From (iii) 



i'(M)=0W„. (25) 



Now, suppose that ip is an orthonormal wavelet. For j G Z let VU,- be 
the closure in L^(R) of the span of {V'j,fc = D^T'^ip : k G Z}. Since is an 
orthonormal wavelet, ((25)) is satisfied. Moreover, the sequence of subspaces 
{V„ : n G Z} defined by ([M)) satisfies properties (i)-(iv). Thus, {V„ : n G Z} 
will generate an MRA if there exists a function (p G L'^ (M) such that the system 
{T''(p : fc G Z} is an orthonormal basis for Vq. In this case we say that the 
wavelet is associated with and MRA or, more simply, that tp is an MRA 
wavelet. 

When ip is an MRA wavelet of L^(]R) with scaling function if the rigid 
functions C of Theorem [TUl and S of Theorem [T51 coincide. Therefore: 

Theorem 20 Letip be an MRA wavelet o/L^(R) with scaling function ip. Then 
there exist rigid functions A, R on J-'L^(R) and B , S on QL'^{M.) such that: 

(i) the initial subspaces of A, B, R are one- dimensional, and S{lo) is unitary 
for a.e. lo G SB; 

(ii) one has 

Wo = J'-^iuj^Ai: k€l) = g-^{S€) , 
Vo=T-^{w''Ri:kE'L) =g-^SH-{dW,P{J)®l^{M: 

(ui) {ip) = T-^{ACj = g-^{i3ej and (ip) = T-^{Ri). 

The MRA wavelet ip and scaling function ip uniquely determine {A, B, R, S) to 
within constant partially isometric factors on the right. 
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Conversely, assume that there exist rigid functions A, R on J-L'^{M.) and B, 
S on QL'^(M.) satisfying (i) and 

T-\u;''Ri: ke1)= g-^{SH-{dB; P(J) © , 

Then the set (A, B, R, S) has associated a unique MR A wavelet ip £ 
with scaling function ip given by 

^ := T~^[Au] = g~^[Bv], if -.^ T~'^[Rw], (26) 

where u, v and w are (constant) normalized vectors in the initial spaces of A, 
B and R, respectively. In such case, (ii) and (Hi) are satisfied. 

Additionally, there exist measurable range functions J and K on ^L^(R) 
such that J{lo), K{lij) are one- dimensional for a.e. lu G 5© and 

Wo = T-'Ttj , Vo = T-^Tlji . (27) 
5.1 Two-scale relations and mirror filters 

Whenever an MRA of X^(R) is given there exists an orthonornial wavelet as- 
sociated with it. Moreover, the wavelet can be constructed explicitly from the 
MRA. This is carry out thanks to the two-scale relations and a pair of discrete 
quadrature mirror filters appears. As we shall see in this Section, these filters 
act just on the spectral model ^L^(M) given in Proposition O 

Let be an MRA of X^(R) with scaling function (p. Since (p G Vq, 

Vo C Vi and {fi^k : fc G Z} is an orthonormal basis of Vi, there exists a 
sequence {hk} G ^^(Z) such that 

fix) =^hk(pi,k{x) =^^hk(p{2x - k) . (28) 

fcGZ fcez 

Equation (j28p is usually known as the two-scale relation of the scaling function 
if. A simple change of variable y = 2^x generalizes this relation to any scale: 

fj.o = ^hk (pj+i,k , (j e Z) . (29) 

fcez 

In terms of the dilation and translation operators, D and T, since ipj,k — D^T^tp, 
relation (|29p is written as 

D^p, = 51 i^^+^TV , (j e . (30) 

Taking j = — 1 in ([50]) and going to the spectral model given in ([3]) , by ([3]) one 
obtains 

[J'D-^p] [lo) = hkLu'' [Tp] {uj) = h{uj) [Tp] [uj) , (31) 
fcez 
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where the two-scale symbol for ip, defined by 

h{uj) i^^/ifetj'^ <EL^{dB), 
fcez 

appears in a natural way. 

In a similar way, if is a wavelet associated to the MRA {^nl^g^j since 
tpj,o G and Vj+i = Vj ® Wj for j G Z, two-scale relations for tp are also 
satisfied: there is a sequence {gfc} G P(Z) such that 

V'j.o = X! 5'=^ 'P]+i,k , U e Z) . (32) 

fc(EZ 

In particular, for j = —1, 

=Y.9kU^' [Tip](cj) = 9{lo) [^^](u;) , (33) 

fcez 

where g is the two- scale symbol for ip: 

The following result is easily proved using the spectral model for T of Propo- 
sition[21 As usual, the symbol _L means "orthogonal to" , the symbol "orthog- 
onal sum" and the overline "adherence of or "closure of. 

Proposition 21 Let {^n}^^^ be an MRA of i^(R) with scaling function ip 
and wavelet ^. Let h and g be the two-scale symbols for ip and ip, respectively. 
Then h and g are elements of L^ satisfying the following properties: 

(a) ± for all n G Z\{0}; 

(b) Lo^'^^giuj) ± g{uj), for all m G Z\{0}; 

(c) uj^"h{uj) ± w2™5(w), for all n, m G Z; 

(d) ®„ez;(w2n/j(^)) 0Q^g2(^2m^(^)^ ^ L^ {OB) . 

Proof: The unitarity of D and the relation T'^D = DT'^'' for all fc G Z give 

This relation together with ([3T|l and (jl]) lead to 

TV-i = ®k&j^{uj-^^ h{uj)[T^]{Lu)) . 
Thus (the symbol means orthogonal complement), 

= .FVo0.FV-i 
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Now, from the two-scale relation ([5^ for ip with j = — 1, 



fcez 

and 

TW-i - efcez(w2feg(^) . (35) 

Since i? is a rigid function and u is a normalized vector of its initial subspace, 
||[^(/?](w)||;2(z) = 1 for a.e uo G 9D and, therefore, equating the right hand sides 
of ([34]) and (|35ll . one gets the result. □ 



A straightforward consequence of Proposition [5T] is: 

Corollary 22 The two-scale symbol h{uj) := X^fcez ^fe^'" satisfies the following 
two equivalent properties: 



hk hk+2n — Sn , (36) 

feez 

|/i(w)p + |/i(-w)p = 2, for a.e. UJ e dB . (37) 



Proof: X^fcez l^fcP = 1 is a consequence of the two-scale relation (|28p and 
the normalization of </5j,fe. That for n G ^\{0} one has '^f^^^hkhk+2n — 
follows from property (a) in Proposition [2TJ The equivalence of ([36|l and (|37p 
follows from writing out the explicit Fourier series for |/i(tj)p -I- a;)p -see 
[H page 137] for details-. □ 

Equivalent conditions to (a)-(d) in Proposition [21] on the functions h and g 
are given by van Eijndhoven and Oonincx in 16 : 

Proposition 23 [van Eijndhoven and Oonincx] Two functions h and g of L'^ (dB) 
satisfy the conditions (a)-(d) in Provosition \21\ if and only if 



(1) g{uj) h{Lo) + g{-Lo) h{-uj) = 0, for a.e. lu e dB; 

(2) The matrix ( !^,^^\ '^''^^ I is invertible for a.e. uj G dD. 

\h{-u) g{~^)J 

Remark 24 From conditions (1) and (2) in Proposition [23] it is possible to find 
g once h is known (or to find h once g is known). For example, it can be verified 
in a straightforward way that possible choices of g are given by 



(w) = uj'^'+'hi-iu) , (38) 



for any m G Z, so that g(w) — X]feez(^l) ^2m+i-fe ^ or 



feez 



gk = i-ir" h2rn+i-k , (keZ) 
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By (|37|) . for the choices of g given in ([38| the matrix ^h{^^^) g[^uj)^ 

unitary for a.e. w G 9D. Thus, conditions (1) and (2) in Proposition [251 together 
with (j38p lead to a pair (ft,, (7) of transfer functions for two (discrete) quadrature 
mirror filters. See [3] Section 5,6] and [8, Chapter 3] for details. 

Remark 25 Property (c) in Proposition [5T] together with ([55]) lead to 

^(-l)i-''^/ife/ii_2„-fc = 0, (nGZ). 

fcGZ 

Remark 26 T/ie iwo discrete quadrature mirror filters {/ifc} and {gk\ act on 
the spectral model J-L'^{W) given in Proposition\^ Indeed, if J-'(p — (p = {(pf"^ }, 
^V^-i.o = ¥5-1,0 = {VP-\.^\f^} and = ■</'_i,o = {['0-i,o]f'-'}, since 

^(^o,fc = - c^'^ <^(c^) - {(^i""'^} , (fc G Z) , 

the two-scale relations ([2^ and (|32p for 7 = — 1 are equivalent to 

Yp^rs^ =Y^h,^-'^ ^ [V^-^r =5:.9.^r^\ (zGl,nGZ). 

On the other hand, relations (pij) and ([55]) can be written as 

(jp_x o(t^) = hiiS) ■ (jpiuj) , ^{yj) ~ g{uj) ■ <f{uj) , for a.e. uj G . 

In general, for j G Z, filtering by {hk} and {(/a:} the coordinates in J^L^(K) of 
(^j,o results in the coordinates of (^Sj-i.o and tpj^i^o, and using the transfer func- 
tions h and g as multipliers on q results in and o; respectively. 

Remark 27 Let {V„}^gg be an MRA of L^(IR) with scahng function ip G Vq 
and orthonormal wavelet ip G Wo = Vi G Vo. Let := 7f+(9D;C) denote 
the space of scalar-valued Hardy functions. Jorgensen [10] considers a spectral 
transform for D^^ of the form 

so that Q^,Vi = i?+(9D, Wo) and Q^,Wq = £ (the subspace of constant func- 
tions), and proves that there exist functions a, bk, b^, {k G N), in such 
that 

ip = a-ip, iff^j, = bk ■ , V>o,-k = K ■ -4^0,-1 ^ (keN), 

where now f :— Q^,f for / G L^(K). Moreover, there is an inner function c, i.e., 
c G and |c(a;)| = 1 for a.e. oj G 9B, such that 

LO c{uj) bk{uj) = blicu) , (w G aD, fc G N) . 

The inner function c serves as a spectral invariant. It separates distinct examples 
and is trivial precisely for the Haar wavelet. 
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Taking into account the two spectral models QL'^{M.) and ^L^(R) of Propo- 
sitions [T] and [2l and the change of representation matrix (a^'n™) , one reaches 
the following identities: 

Proposition 28 Let {V„}^g^ be an MR A of L^(R) with scaling function ip 

and = (p = {i^l"-*}. If h{uj) — X^feez^fc'^'' ^■^ ^^^^ two-scale symbol for if, 
then, for every pel and q ^1, 



s,j,m i^n k {^C\\ 

E\ ^ s.j.m \ ^ /t s. 7,771 — 1 , 7 s.j.m — 1\ '-('".) ^ 

k s,j,m i,n 

Moreover, the function ip G L^(R) defined by Tij: = {V'!"''}, where, for every 
p G I and q & "L, 



^i'' = E E<'r"' E(-i)'~''^i-^^i""'^ 

sj,m A; 

E\ ^ s.j.m \ /i s.7,m— 1 7 s.7,m— 1\ ^(^^) ^ ' 

s,j,rn i,n 

is an orthonormal wavelet associated with the MR A |Vr>| 

Proof: The two-scale relation (p-i.o — J2k Vo.k on TL'^(R) looks like 

,p ^ Dih ■ >p) = D(^J2hkf'' = Y,hkDf''g>, (41) 

fc k 

where £) := TDT'^ and f := TTT-^. By dH), the first hne in JSll) is the 
counterpart of (|1T|) for the coordinates {(^|"^}. Moreover, since DT"^^ =T^ t) 
for every k £ Z,, 



(p^Y.khkDf^(p =Ek [h2k D f^'' + h2k+i D T^'-^+i] (p 
= Ek T'^D [h2k + h2k+i f] (p . 



(42) 



The second line in follows from dH]), (O and ©. Finally, that the function 
given by P0)l is an orthonormal wavelet associated with the MRA {Vn}^^^ 
is a straightforward consequence of Remark [24] with m = 0. □ 

Remark 29 Not every pair {h,,g) of functions of l?(dV) verifying conditions 
(a)-(d) of Proposition[2T]or conditions (1) and (2) of Propositionl^ comes from 
an MRA. Under additional conditions, Cohen [2] and Lawton [9j gave equivalent 
conditions for it. This question together a deep analysis of relations (p9)) and 
(PD|) shall be considered elsewhere. 
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A Some bases for the spectral models 

In order to obtain the respective spectral models given in Propositions [1] and 
[2] for the dilation and translation operators, D and T, one must consider or- 
thonormal bases (ONB) of L'^i'R), 

{-ftri™' : j e J, m e Z} and : i G I, m G Z} , 

such that {K^f^j ■■ j el} is an ONB of L^[±l, ±2), {if ^ : i e 1} is an ONB of 
[0,1) and 

The change of representation between the two models is governed by the matrix 
{al'n"^) whose elements satisfy 



A.l Exponential bases 

For the exponential bases 



K^±lix) e^"^'^ and L^°\x) - -'^^^'^^ 



e 



one has I = J = Z and the change of representation matrix (a^'-Jj™) is given by: 

-^TT-, ^^2^^k2-"^ if s = +, m > 0, /c - j2" ^ 0, 

_ 27r(fc-j2™) ^ ' : .V 7- , 

-fc,o ^ 2-™/2, if s = +, m > 0, k- j2" = 0, 

0, otherwise. 



• 2m/2 -27ri/c2"™ 

— — -— (e"2'^*'^-2 " - 1), if s = m > 0, fc - j2" ^ 0, 

27r(fc-j2") ^ ' ^ J f , 

2"™/^ if s m > 0, fc - j2'" = 0, 

0, otherwise. 

s,].m ^1 1, if s = +, TO = 0, A: = j, 
'^'^ 10, otherwise. 
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1, if s 



TO = 0, k = j, 



'^'"^ 1 0, otherwise. 
For p = 1, 2, . . . and g = 0, 1, . . . , 2P - 1, 



^ 2P/2 ^-2-Kij2-Pq 

27r(A:2P - j) 
2-f/2, 
0, 



-27rij2-»' 



1), if s = +, m = -p, A;2P - i 7^ 0, 

if s = +, TO = -p, A;2P - i = 0, 
otherwise. 



at'"' 



fe,-2P-g-l 



= < 



27r(A;2P - j) 
2-f/2, 
I 0, 



„27rij2-P 



1), if 



-, TO = -p, fc2P - i ^ 0, 



if s = — , TO = — p, k2''' — j — 0, 
otherwise. 



A. 2 Haar bases 

Let tp be the Haar wavelet and (p the corresponding scaUng function, 

= X[o,i/2) - X[l/2,1) > "/^ = X[0,1) ' 
where the characteristic function x is defined by 



1, if a; e [a, b), 
0, otherwise. 



As usual, for each j, fc G Z, let ipj^k ■= D^T'^tp and (/?j,fe := D^T'^ip. 
Consider the Haar bases defined by 

Kf^2P+q ■= i for p = 0, 1, . . . ; g = 0, 1, . . . , 2^ - 1 ; 

■f'^-.o := ; 

^i^i"2P+, := V'p -2f+i+g , for p = 0, 1, . . . ; g = 0, 1, . . . , 2P - 1 : 

^0 '^0,0 I 



(0) 
^2P+(j 



:= Vp,,, forp = 0,l, g = 0,1,..., 2^-1. 



Then I = J = N U {0} and the change of representation matrix (a*'^''") is as 
follows: 



For n = 0, i = 0, 



0,0 



2""/^ ifs = +, j = 0, TO>0, 

0, otherwise. 
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For 71 = 0, r = 0,1,2,..., 

r -2-1/2^ if s = +, j = 0, m = r + l, 
a^;^;™ = <^ 2('-™)/2^ if s = +, j = 0, m > r + 1, 
[ 0, otherwise. 

For n = 0, r = 0, 1,2, . . . and t = 2^ + 9 (with < p < r and q = 0, 1, . . . ,2^-1), 

^s,j,m ^ / 1' if s = +, j = t, m = r - p, 
2'-+t,o 1 Q otherwise. 



For n = 1 and every « > 0, 



1, if s = +, j = z, TO = 0, 
0, otherwise. 



For n > 1 and i > 0, with n = 2'' + t; (u = 1, 2, . . ., w = 0, 1, . . . , 2" - 1) and 
i = 2'- + t (r = 0, 1, . . ., t = 0, 1, . . . , 2'' - 1), 

^ / 1' ifs = +, j = 2'-(2" + i;)+t, TO=-u, 
2'-+t,2"+u ~ 1^ 0, otherwise. 

For n > 1 and i = 0, with n = 2" + i; (u = 1, 2, . . ., w = 0, 1, . . . , 2" - 1), 



an 



2""/^ if .S = +, j = 0, TO = -u, 

_ J 2(P-")/2, if s = +, j = 2P + K2"-P] 

for < p < w, TO = — M, 
0, otherwise. 



where [•] denotes "entire part of an 

V - 2"-P[i;/2''-P] 



w{u, v,p) 
For n = —1, i = 0, 



2^— p— 1 



(43) 



2-"/^ ifs = -, j = 0, m>0, 
1 0, otherwise. 



For 71 = -1, r = 0, 1,2, , 



2"^/^ if s = -, j = 0, TO = r + 1, 

1,-1 = ^ -2(''-")/2, if s = -, J = 0, TO > r + 1, 
0, otherwise. 



s.i.m 

"2 



For 71 = -1, 7' = 1, 2, . . ., < p < r and g = 0, 1, . . . , 2P - 1, 



2.+i_2P+i+g.-i ] 0, Otherwise 



1, if s = — , j = 2P + 771 = r — p, 



^Taking into account the binary expression v = Xlfe^o^fe^'', with tj, = or 1, one has 
«;{n,t;,p) = i„_p_i and [v/2^-P] = Eft=Lp *fc 2'=-("-P) . 
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For n — —2 and every i > 0, 

1, if s = j = i, m = 0, 
0, otherwise. 

For n < -2 and z > 0, with n = -2"+^ +v (u 1, 2, . . ., w = 0, 1, . . . , 2" - 1) 
and i = 2"^ + i (r = 0, 1, . . ., t = 0, 1, . . . , 2'' - 1), 

s,j,m ^ / 1, if s = -, j = 2'^(2" + i;)+t, m = -w, 

"a^+t -2-+i+« I 0, otherwise. 

For n < -2 and i = 0, with n = -2"+^ + v (u = 1, 2, . . ., w = 0, 1, . . . , 2" - 1), 

2"''/^ if s = j = 0, m = -u, 



for < p < u, m = — M, 
0, otherwise, 



where 'w{u,v,p) is given by ((43|) . 
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